will be published in Journal of Applied Physics 

Wave excitations of drifting two-dimensional electron gas under 

strong inelastic scattering 

V. V. Korotyeyev, V. A. Kochelap 
Department of Theoretical Physics, Institute for Semiconductor Physics, 
Pr. Nauki 4I, Kiev 03028, Ukraine 

L. Varani 

Institut d'Electronique du Sud, CNRS UMR 5214, University Montpellier 2 France 

Abstract 

We have analyzed low-temperature behavior of two-dimensional electron gas in polar heterostruc- 
tures subjected to a high electric field. When the optical phonon emission is the fastest relaxation 
process, we have found existence of collective wave-like excitations of the electrons. These wave-like 
excitations are periodic in time oscillations of the electrons in both real and momentum spaces. 
The excitation spectra are of multi-branch character with considerable spatial dispersion. There 
are one acoustic-type and a number of optical-type branches of the spectra. Their small damping 
is caused by quasi-elastic scattering of the electrons and formation of relevant space charge. Also 
there exist waves with zero frequency and finite spatial periods - the standing waves. The found 
excitations of the electron gas can be interpreted as synchronous in time and real space mani- 
festation of well-known optical-phonon-transient-time-resonance. Estimates of parameters of the 
excitations for two polar heterostructures, GaN/AlGaN and ZnO/MgZnO, have shown that exci- 
tation frequencies are in THz-frequency range, while standing wave periods are in sub-micrometer 
region. 

PACS numbers: 72.30.+q, 72.20.Ht, 73.63.Hs 
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I. INTRODUCTION 



In semiconductor materials and heterostructures at low temperature, when ^-'^^I^bT ^ ^ 
(o'o, and T are the optical phonon frequency, the Boltzmann constant and the temper- 
ature, respectively) the absorption of optical phonons by electrons is practically absent. 
While optical phonon emission can be the dominant scattering mechanism for hot electrons 
if electron-optical phonon coupling is strong enough. In samples with high electron mobility, 
where the electrons undergo only weak quasi-elastic scattering at a low field, the dynamics 
of an electron subjected to a steady state high electric field Fq is the following. The elec- 
tron is almost ballistically accelerated by the field until reaching the optical phonon energy. 
Then, the optical phonon emission occurs so that the electron loses practically all its energy 
and stops. The processes are repeated in time. This dynamics gives rise to temporal and 
spatial modulation of the electron momentum, p, and velocity v with characteristic time 
period, tf — pq/cFq, and length, Ip = eFQTp/2m* = hcuo/eFo, where po — \/2m*l%u}Q, e is 
the elementary charge and m* is the electron effective mass. 

For this kind of cyclic motion, the electron momentum space can be divided into two 
regions. The first is the passive region, where the electron energy £'(p) < fiuiQ. In this 
region inelastic scattering by optical phonons is almost absent, so that electron scattering 
is practically elastic and can be characterized by an elastic scattering time, Tp. The second 
region is the active region, E{p) > hcuo, where optical phonon emission is the dominant 
process, i.e. 

(1) 

with being the optical phonon emission time. Pronounced cyclic motion is possible, 
when 

C«T^<T^, (2) 

i.e., the time of flight of the electrons through the passive region is shorter than the scattering 
time, while penetration of the electrons into the active region, p — po = Ap, is sufficiently 
small, 

Ap~po-^<Po- (3) 

Actually, the latter inequahties deflne a range of the electric flelds, where this kind of periodic 
motion can occur. The critical electric field necessary for the onset of the cyclic electron 
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motion is estimated to be 



^cr-—. (4) 

eTp 

The presented scenario is essentially single electron physical picture valid at low or modest 
electron concentrations, when e-e collisions do not destroy the cyclic motion. 

The possibility of cyclic acceleration-stop electron dynamics due to strong inelastic scat- 
tering by optical phonons was predicted six decade ago by Shockley.- Experimental evidences 
of the cyclic dynamics in real space were found by analyzing low temperature I-V characteris- 
tics of short diodes made from different polar materials: InSb,- InGaAs,- GaAs,- and InP.- 
Tens of cycles were identified at low temperatures. Very recently, comprehensive Monte- 
Carlo simulations of electrically biased short InN and GaN diodes have demonstrated the 
formation of stationary one-dimensional gratings of electron concentration and velocity at 
nitrogen temperature.-i^ 

In time (frequency) domain, the cyclic dynamics gives rise to a resonance phenomenon at 
the transit-time frequency Up = 2tt/tf, frequently called optical phonon transient time res- 
onance (OPTTR). The OPTTR induces a number of interesting effects (for review see Refs. 

nci 

8|j9|). The most important manifestation of this phenomenon is the possibility of microwave 
amplification and generation in the sub-THz and THz frequency regions. The microwave 
generation based on OPTTR phenomenon was observed experimentally in InP samples for 
the frequency range 50 to 300 GHz.^^ For other indicated above bulk-like materials, the 
OPTTR generation was studied theoretically in details.-i^ 

Discussed cyclic electron dynamics should be even more pronounced for low-dimensional 
electron systems. Indeed, for the latter the onset of the optical phonon emission is much 
sharper because of specific features of the density of low- dimensional states. For comparison. 



the optical phonon emission rate I/t^^ in bulk crystals is proportional to ^yE{p) — Tiujq (0 < 
i?(p) — TujJq <^ huo), while for two-dimensional carriers, l/r^^ oc TilE^p) — huo], where H[x] 
is the Heaviside step-function. In addition to a sharper threshold of the optical phonon 
emission, high mobilities and, thus, quasiballistic motion at energies E{p) < Huq can be 
easily achieved for low- dimensional electrons, while their concentration can be controlled to 
avoid e-e scattering. These important qualitative conclusions are supported by papers,—"— 
where, by the example of GaN quantum well heterostructures, it was shown the preference 
of their use to reach the cyclic electron motion. Note that effects similar to OPTTR are 
predicted for novel one- dimensional nanostructures - carbon nanotubes.— 



Summarizing this short review, one can conclude that the cychc electron dynamics is 
a quite general low-temperature phenomenon characteristic of many polar materials and 
heterostructures subjected to an electric field. It is believed that the phenomenon may 
lead to different effects, the most representative of them are spatial grating of electron 
concentration and velocity in finite size diodes and OPTTR resonance in the time/frequency 
domain for bulk samples. 

In this paper, by investigating the case of two-dimensional electrons, we show that under 
the conditions of the cyclic electron motion, a novel type of weakly damped excitations of 
drifting electron gas exists. These excitations, quite different from well known plasmons, are 
periodic in time and in real space oscillations (waves) of electron concentration and charge 
both synchronized with electron redistributions in the momentum space. Their frequency- 
wavevector relation has infinite number of continuous branches, 00^ {q), with q being the wave 
vector of excitations and k = 0, ±1, ±2.... The excitation damping is especially weak or even 
absent, when the frequency and/or the wavevector are multiples of up and/or qp = 2tt/Ip, 
respectively. Particularly, for q multiple of qp standing waves {u =0) without any damping 
appear. 



II. MODEL AND EQUATIONS 

It is well established that, for conditions when the cyclic electron motion can occur, steady 
state and high frequency electron transport can be described by semiclassical Boltzmann 
transport equation (BTE) for the electron distribution function in the momentum space. 
Under conditions ([2]), the electron distribution is mainly grouped in the passive region and 
essentially elongated along the electric field (the so-called streaming distribution). The BTE 
can be solved numerically by Monte-Carlo simulations or analyzed analytically with the use 
of same approximations adequate to the physical picture. These approximations include: 
the Baraff approximation^ to take into account the anisotropic part and setting to zero 
the isotropic part of the distribution function in the active region because of the short 



emission time r^™ (see discussions in Refs. 
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171 for bulk materials, and in Refs. 



11 



for 



two-dimensional systems). In this paper we will use the analytical approach. 

Consider a uniform two-dimensional electron gas confined in narrow quantum well layer 
at 2; = (the quantum well thickness is negligible in comparison with the characteristic 
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lengths under consideration). Introduce the reference frame with X,Y-a.xes in the plane of 
the quantum well layer and Z-axe perpendicular to this plane. The applied electric field 
is along the X-axe. Nonequilibrium electrons are described by the distribution function 
$ = ^{px,Py,r,t), which is, generally, dependent on the momentum, p = {Px,Py}, the 
coordinate, r = {x,y}, and time, t. Let be the areal electron concentration given by 
doping. Then we define the distribution function in such a way, that A^<l'(p, r, t)dpdr is the 
number of electrons in the phase-space elementary volume dpdr located at the phase-space 
point {p,r}. The BTE for $ is 

^ + - = Leim + Lo,m , (5) 

at m* or op 

where the collision integrals Lez{^} and I/op{$} describe the quasielastic scattering and the 
inelastic scattering by optical phonons, respectively. 

An excitation of the drifting nonequilibrium electron gas can be considered as a per- 
turbation of the stationary distribution with appearance of an additional field. Thus, we 
present the distribution function and the field as stationary and perturbed contributions: 
$ = $0 + $ and F = Fo + F, where F = {-Fq, 0, 0} and the field F is defined by the Poisson 
equation: 

ineN f 

divY = -^^5[z\ / rfp$, (6) 



where Kq is the dielectric constant of the quantum well surroundings; 5[z\ is the Dirac 
function. The perturbations are assumed to be dependent on time t and x coordinate as 
follows: $(p, a;,t) = (j){p)exp[iqx — iut], Y{x,z,t) = i{z) exp[iqx — iut]. In this case the 
field induced by the perturbed electron gas has two components: f{z) = {fx{z),0,fz{z)}. 
The Fourier components 0(p), fx{z),fz{z) are dependent on u and q. 

The BTE in its general form of Eq. ([5]) can be simplified by adapting it to the problem 
under consideration as follows. Both contributions to the distribution function, $o and 0, 
are approximated as sums of isotropic and anisotropic parts:— 

where the two-dimensional vector p is presented by its modulus p and the angle 6 to the 
X-axe. 

The next simplifications are the following. In the passive region, for elastic momentum 
scattering it is easy to obtain Lei{(po} = l/ei{0o} = 0. While in equations for anisotropic 



parts, 4>i,4>i, the elastic scattering integrals produce the relaxation terms — {(/)i, </)i}/rp. We 
suppose that the elastic scattering is due to interaction with acoustic phonons. Then, for two- 
dimensional electrons, Tp does not depend on the momentum/energy. In the active region, 
the collision operator describing optical phonon emission can be estimated as Lop^ ~ ^/t^^, 
where the optical phonon emission time is the shortest time in the system (r^™ — > 0). 
Together with inequality ([1]), this provides the following conditions^ for the isotropic parts 
of the distribution function: 0o(p > Po),(t^o{p > Po) ~ 0. In contrast, the anisotropic parts 
of the distribution functions, 0i, 0i, are finite at p = po due to field induced electron stream 
in the momentum space. However, as discussed above, the number of streaming electrons in 
the active region is rapidly decreasing because of the optical phonon emission. As a result 
of these simplifications, the BTE is reduced to a system of coupled ordinary differential 
equations for 0o, 0i, 0o, 0i in the interval < p < pq. 
For the steady-state problem these equations are: 



^dcpo I 1 d{pct>i) 
dp p dp 







(8) 



Two first-order differential Eqs. (|8]) have to be supplied by two conditions to determine 
integration constants. One of these can be obtained from the normalization of the steady 
state function. We set / ^o{p)dp = 1, then we have 



p dp [27r0o + 



The second condition follows from the above discussion: 



PO 



(9) 



00 bo) = 0. 

The spatial and temporal problem is described by equations: 



i 

+ 



-a;(27r0o + 0i) + 5<^i 



+ e 



-tu 



Idipjnl 
p dp 



Fq d(p(t>i) 
p dp 

feFn 







dp 



(10) 



(11) 



In order to obtain the "boundary" conditions for the latter system, we use the continuity 
equation, that follows from the BTE in its general form of Eq. ([5]): j dp + = 0. 

Using our notation, it reads 



■po 



p I 



pdp [-w(27r0o + 0i) + 0i] = . 
m* 



(12) 



Then with the use of the first equation from (fTTi) . we obtain: 

pi>i\p^o - PoMpo) = . (13) 

The second boundary condition is similar to Eq. ( ITOj) : 

0oM = O. (14) 

The ac-electric field component, fx, that enters into the right hand side of the second 
equation of the system fllip . is defined via the Fourier component of the electrostatic poten- 
tial, ip: fx = —iq(p. The latter obeys the Poisson equation: 

- gV = / dpp 27r0o + 01 15 

with the boundary conditions (p{z — )■ ±oo) = 0. These relations finalize the mathematical 
formulation of the steady-state and wave-like excitation problems for the streaming two- 
dimensional electron gas. 



III. SOLUTIONS OF EQUATIONS 

To proceed further it is convenient to introduce the following dimensionless variables: 

X = ^,T = —,P = ^,n = urp,T = ^,Q = qlp. (16) 
If tf Po Tp 

We substitute Pq0o ^ 0o ) Po4>o ~^ 4>0y etc. Thus, the functions (po, 0i, 0Oi 4>i ^^^^ ^ow 
dimensionless. 



A. Solutions for steady state problem 



The solutions of Eqs. (l8l)- (]T0l) for the steady state problem are 

— 2r 

MP) - pl^i (18) 

Eqs. ( IT71) describe the electron distribution quite well almost in the whole momentum 
space, except a narrow interval corresponding to the electron penetration into the active 
region given by Eq. ([3]) and the respective interval nearby p = 0.— 



Having /o and /i one can calculate the average electron energy and the drift velocity: 

„2 8 -I- 

^2m* 12[2 + r] ' ^ ^ 



7 ^m* [2 + r] m* ' ^ ^ 



For the case of the strong inequality ([2]), we obtain the characteristics of the "ideal" electron 
streaming regime: = ^c^o/S, fdr = Po/2m*. 



B. Solutions for the space and time dependent problem 



Using the dimensionless variables of Eqs. f|T6|) one can rewrite equations Eqs. ffTTj) in the 
form: 

d(P4>i) 



I 
i 

+ 



n 27rP0o + P01 + 2QP 



+ 



dP 



0, 



[-n - 2F) P01 + 2gP 7rP0o + P01 



TT- 



d{P4>i) 
dP 



(21) 



" Fo dP ■ 



The following substitutions 



7rP0o = Xoexp [2 (^]P-gp2)] , 
P0l = Xlexp[^(^]P-gp2)] , 



(22) 



transform Eqs fl2Tl) to the simpler form 



% = 2znxo 



X exp [-Z (1]P - gp2)] . 



(23) 



The latter system of two equations can be reduced to a single equation for the function, 
Xi(P), that we present as 



(fxi 



+ 



b- 



1' 
P 



dx 



dP 



l + cxi = ^R{P\n,Q) 



Pr 



(24) 



where we designate 



b = 3iQ, c = 2iQr, 
R{P \Q, Q) = 27TiQP^ exp [-i {VlP - QP^)] . 



(25) 
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Note that, in the left hand side of Eq. (^^, the parameters, b, c, are dependent only on fl, 



while the right hand side parametrically depends on both Q and Q. From Eqs. ([13 
and the first equation of the system (125]) it follows the boundary conditions to Eq. 

xm = Xi{l)exp[t{Q-Q)] , 
S^(1) = 0. 

The general solution to Eq. can be presented as 



(26) 



Xl(P) = Ci^l(P)+C2^2(P) + 



fx 

To 



dP'-^>2{P) j 



Wn{P') 



dP' 



(27) 



Wn{P') 

where \l/i(P), \1'2(P) are two independent solutions of the congruent homogeneous equation, 
Wn{P) is the Wronskian of this equation, Ci, C2 are arbitrary constants. The functions 
\E'i(P) and '^2{P) can be expressed via the Kummer functions of first and second kinds,— 
M{u,w,z) and U{u,w,z), respectively: 

^i(P) = p2exp(-aP)M[/3,3,eP] , 
^2(P) = P' exp{-aP)U[f3, 3, eP] , 

with e = Vb"^ - 4:C , a = {b + e)/2 and (3 = 3/2- b/2e. 

For the problem under consideration, the constants Ci and C2 should be found by using 
the boundary conditions ( 126|1 . The latter result in the following equations: 



Ci {^i(O) - ^i(l) exp [i {n -Q)]} + C2 {^2(0) - ^2(1) exp [i {n - Q)]} 

fx 



[^l(0)J2(0)-^2(0)Jl 



(28) 



(29) 



-1 RiP')^l,2(P') 



where Ji,2(P) = Jp ^^(p,) 
side of Eq. ( l29l) , is defined by the Poisson Eq. ( !T5l) : 



dP'. The ac electric field, fx, that appears in the right hand 



— = sgn{Q)—— 



dPP 



sgn{Q) 



2meN 
1^0 Po 



dP 



Xi 



i dxi 

n dP 



exp [i {QP - QP^)] 



(30) 



where sgn{Q) stands for sign of Q. Substituting f lTH) to Eq. fl5U]) one can find fx as a linear 
homogeneous form of the constants Ci, C2: 



^=U{A,C^ + A2C2) 



(31) 



where Ai^2 = ^i,2{^ij ^2!^, Q} are expressed via bulky integrals from some combinations 
of the functions and "^2- We introduce the parameter 

Af=^, N^H^^, (32) 
Nch ' ''^ 27re ' ^ ^ 

which has the meaning of a dimensionless electron concentration. From the structure of 
Eq. fl29l) it follows that this parameter determines the coupling of the dynamics of individ- 
ual electrons and their collective motion arising due to the self-consistent electric field, f. 
Eqs. ( l28l) . ( 129|) supplemented by the relationship ( 1311) compose the system of homogeneous 
linear algebraic equations for the constants Ci, C2. The solvability condition of this system 
is determined by the zero of the associated determinant, A{Q, Q) = 0. The latter condition 
gives the dispersion relation, fl{Q), for required excitations of the drifted electron gas. 

C. Frequency and wavevector dispersion of electron conductivity and dielectric 
permittivity of the electrons. 

The solvability of the obtained equations can be also presented in a more usual form. 
Indeed, the high frequency current can be calculated as 

i = -^ fdPP%{P), (33) 

III' Jo 

where parametrically dependent on Q and Q function 0i(-P) is the solution of Eqs. (!22|) . 
( l27ll . ( !28l) . ( l30l) obtained at a given fx/Fo- Then, the high frequency conductivity is 

a{n,Q) = t= {n, Q) ^ ei^oN^u {n, Q) , 

fx 'f^^O i'O 



i,{Q,Q)^-^ f'dPP^MP), (34) 

Jx Jo 



where /xq = eTp/m is the low field mobility, z/ (17, Q) does not depend on /x/Fq and has 
a meaning of dimensionless electron mobility dependent on both the frequency and the 
wavevector. 

With the help of Eq. ( 134|1 . one can calculate the dielectric permittivity of the electrons: 

e (fi, g) = 1 + ^2ArM ^ (fi, g) . (35) 

It can be easily proved that the usual requirement 

e{n,Q)=0, (36) 
10 



that defines the possible collective excitation modes of the electron gas, is identical to the 
above obtained solvability condition. 

For what follows, it is useful to discuss briefly some properties of the high frequency 
mobility, z/(i7,g). Fig. [T] shows the results of calculations of z/(fi, Q) for F^/Fcr = 3. The 
OPTTR becomes apparent near Q = ±2tt. Indeed, one can see typical frequency dispersion 
of the conductivity in narrow frequency intervals and, particularly, a change of the sign of 
the dynamic conductivity (i?e[z/] < 0). The latter feature has been in focus of previous 
studies.-i^i^ii-^ At Q 7^ 0, drift of the electrons is revealed: the characteristic feature arises 
at i7 ^ VdrQ, where Vdr is the dimensionless drift velocity defined as Vdr = 2vdrm/po with 
the field dependent Vdr given by Eq. ( EUj) . Besides, a negative dynamic conductivity occurs 
at the Cerenkov region, < < VdrQ, while the OPTTR frequencies are shifted by factor 
VdrQ- Note, the absolute values of z/(n, Q) are relatively small. This is typical for conditions 
of the streaming regime, when the drift velocity saturates and the differential mobility tends 
to zero. 

IV. SOLUTIONS AT ^ 0. 

First, we study the space and time dependent solutions neglecting self-consistent electric 
field f. According to the relationship ( 13T1) . this corresponds to the limiting case 



In this limit, the electron gas excitations can be thought as phased motion of the electrons 
in the real space and the momentum space under the conditions of cyclic dynamics of 
individual electrons. Importantly, all results for this limit depend on the single parameter. 



For such a case, from Eqs. (1280 and (129|1 we immediately obtain the condition of solvability 
of these equations in the form 



Ar< 1. 



(37) 



cr 



Tp/Tp = 1/r. 



Ao (r],Q) = ^i(o) 



{1)-Wn{l)exp [i {n-Q)]=0, 



(38) 



dp 



dp 



that can be rewritten as 



exp {-iQ) 




exp{-in) = s{n',n"\Fo). 



(39) 
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In the latter form, the left hand side depends on Q, while the right hand side is a function of 
Q. We will look for wave-like excitations with a real wavevector, Q, while the frequency can 
be a complex value, Q = Q' + iQ". To find fl' and fl" at a given Q we have two equations 



where k is an integer, k = 0,1,2.... Here we introduce the absolute value, abs[S], and 
the phase, arg [S], of the complex function S{Q',Q"). Eq. f l40|) defines the damping of the 
excitations, Q", as a function of the frequency Q'. Then, Eq. (HT]) defines implicitly the 
dispersion relation of sought-for wave excitations, Q'{Q) for the limiting case of Eq. f l37|l . 

To understand the nature of the solutions of Eq. (HO]) , one should analyze the following 
function of two variables Q{Q',Q") = abs [S{Q' — 1. This function parametrically 
depends on a single value: the dimensionless field Fq/Fct- In Fig. [2], a typical density plot of 
Q{Vl' , Vt") is presented for Fq/ F^r = 3 (see below discussion of actual parameters for particular 
materials). As seen from this figure, there are two kinds of solutions, Vl" = Q" [Q']. The first 
is characterized by small Q" and, particularly, Q"{Q' = 0) = 0. The second kind of solutions 
gives Q" of the order of Q'. Below we will concentrate on the first kind of solutions, because 
the electron gas excitations corresponding to such solutions are weakly damped. 

From Eq. (14T!) . it follows that there is an infinite set of branches of the dispersion relation 
corresponding to different values of the integer k. We will call the branch of /c = as the 
acoustic-like dispersion branch. From Eqs. ( HOj) and fHTj) . for the acoustic branch one can 
find a function Q = Q^{Q). Then, any other 'satellite' branch can be obtained by a simple 
translation in the {fi, Qj-plane along the Q-axe: Q = Q''{Q) = Q^{Q ± 27Tk). The "basic 
vector" of translation is Qo = 2tt. Its dimension value is 2tt/If. In the real space, this value 
corresponds to a wavelength of the excitation strictly equal to the characteristic spatial 
period of the cyclic electron motion, Ip. For the acoustic branch, we find that ~ VdrQ at 
Q — > 0. At finite Q, the acoustic branch can be obtained numerically: the result is shown 
in Fig. [3] (a). In this figure, we presented also two satellite branches with k = —1, 1. The 
latter branches have finite frequencies at Q — ?■ 0, i.e., they can be thought as of optical type. 
Interestingly, the frequencies of the satellite branches at finite wavevectors Q = Qok are 
exactly zero. This means that the corresponding electron gas excitations are standing waves. 
For convenience in Fig. |3] (b), we show also the imaginary part of the frequency, Q"{Q), for 



Q 



abs [s{n', n")] = 1 , 

arg [S{n', fi")] ± 2nk , 



(40) 
(41) 
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each of the branches. The absolute values of Q" are small in comparison with Q'M The 
solutions corresponding to the standing waves have zero damping, Q"{kQo) = 0. 

For the latter particular case {Q = 0, while F 7^ 0), the solutions can be found in an 
explicit form. Indeed, from Eqs. fl25]) we find Xi = C = constant , Xo = —CVPlnP. Then, 
from the boundary conditions (12^ we find the allowed Q = 2TTk ,k^0. According to 
Eqs (I22I) . contributions to the distribution function are 

CT 

00 = InP X exp{-i2'KkP^) , (42) 

TT 

c 

0^ = _ X exp{-i2TikP^) , (43) 

Interestingly, the Fourier component of the perturbed electron density corresponding to 
these solutions is not zero, while the perturbation of the average electron velocity is zero, 
i.e. there are no electron fiuxes in real space. 

Additional analytical results can be obtained for the ultimately simple case adequate to 
F = (absence of any scattering in the passive region). In this limit, from Eqs. f l23|) . we 
obtain xo = 0, Xi = = constant. The corresponding contributions to the distribution 
function are: 

00 = 0, 01 = ^exp[iinP - QP^)] . 

Then, from the boundary conditions (|26l) we find exp[i{Q — Q)] = 1, i.e., the dispersion 
relation gives an infinite set of straight lines: ^Y^oiQ) = Q + 2TTk , /c = 0, ±1, ... (F = 0) . 
At nonzero but small F, this dispersion relation can be corrected by using the perturbation 
method: Q = f^r=o + ^^iQ + 271 k), where the function B{X) is defined as 

Jo Ji y 

Near the points Q = — 27rfc, where f^r=o small, we obtain approximately 

n^{Q + 27ik) ( 1 - - ) - t-{Q + 2nk)^ . 



3 

Thus, electron momentum dissipation in the passive region leads to the damping of the wave 
excitations, as well as to renormalization of the excitation frequency. However, the standing 
waves are not damped. 

The found dispersion relation and solutions can be explained by time-dependent phased 
motion of many electrons in both the real space and the momentum space. The properties 
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of such a phased motion are very distinct for different Q and k. Indeed, for the acoustic 
branch a.t Q, Q the solutions correspond to the simple drift of a perturbed electron 
concentration with velocity Vdr (this type of solutions can be obtained, particularly, in the 
simplest hydrodynamical models). At /c 7^ and small Q (optical-like branches), the time- 
dependent phased motion occurs mainly in the momentum space. At 7^ 0, Q ~ kQo, the 
electron motion in the real space is phased with that in the momentum space, but it is time 
independent. 



V. SOLUTIONS AT TV OF A FINITE VALUE. 

For the case of a finite value of A/", when collective charge effects are involved, we per- 
formed a numerical analysis of the equations. Before presenting this analysis, we note that, 
due to the scaling given by Eqs. (IT^ . the results obtained for A^ = were dependent only on 
the dimensionless parameter Fq/Fct. Now, for a given N we calculate uiVt^Q) and eiVt^Q) 
according to Eqs. then we solve Eq. (15^ for the complex frequency at a given 

real wavevector Q. Representative examples of these calculations are given in Fig. HJ where 
the real and imaginary parts of Vt are shown as functions of Q for three branches of excita- 
tion spectra at F^/F^r = 3 and A/" = 7.4. Comparison of the results for the real part of the 
frequency, Q', with those presented in Fig. [3] shows that, at a finite A/", noticeable changes 
occur for the acoustic branch of the excitation spectra, while changes of the optical-like 
branches {k = ±1) are much smaller. Particularly, there is no effect of A/" on the existence 
of the standing waves {Q = 0), which have the same spatial periods {Q = 2'Kk). A more 
detailed evolution of the acoustic-like branch when M varies is illustrated in Fig. [5] (a). 
Comparing the damping of the excitations, fi", presented in Fig. H] (b) one can make similar 
conclusion: at finite values of M the damping of the acoustic-like excitation branch increases 
significantly, while that of the optical-like branches is less affected. Specifically, the standing 
waves remain strictly undamped. In Fig. [5] (b) the damping of the acoustic-like branch is 
presented for different A^. 

Remarkably, the solutions for the standing waves can be obtained exactly for arbitrary 
N . The solutions include, particularly, the expression for the amplitude of the electric field, 
fxi associated with space and time dependent redistribution of the charged carriers. These 
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solutions are: 



n = 0,Q = 27Tk 



InP X exp{—i2TTkP'^) 



'0 = 



vr 2 + T + AiTAfl2ik) 




>! = — X exp{—i2TckP'^) 

[2 + r] [/i(A;)-2rJ2(A;)] 
2 + T + AiTXl2{k) 



(44) 



(45) 



Here C is an arbitrary constant, functions Ii{k) and /2(^) are 




dPe 



, /2(A:) = / 



Jo 



dPPlnP e 



As indicated above, the found solutions correspond to a time-dependent motion of the 
electrons in real space highly correlated with that in momentum space. Spatial charge ef- 
fects are related to electron behavior in real space. To analyze them, one can introduce the 
concentration of the "isotropic" electrons, Nq, and the concentration of the "streaming" elec- 
trons, A''i, which are determined through the disturbed distributions 0o and 0i, respectively. 
The total perturbed concentration of the excitation wave is = A^^o + A^^i. To understand 
the main characteristics of the time and real space dependences of these concentrations we 
neglect the small damping of the excitations, so that the dependences are 



Here we take into account that a given spectrum branch number, /c, and wavevector, Q, 
determine the excitation frequency, f2'^(Q), and the particular solutions. Then, 5q and 5i 
are phase shifts dependent on Q and k. 

In Fig. [6] the relative amplitudes of concentrations of "isotropic", and "anisotropic" elec- 
trons, and their phase difference 5{Q, k) = 6o{Q, k) — 5i{Q, k) are shown for three excitation 
branches presented in Fig. |H The results illustrate useful symmetry properties of these func- 
tions. Indeed, for the acoustic branch (A; = 0), the relative amplitudes, A'o,i((5, k)/N{Q, k), 
are even functions of Q, while the phase difference is an odd function of Q. For a given 
optical branch, the relative amplitude is a strongly asymmetrical function of Q. However, 
there are simple relationships between the characteristics of two branches with ±k. At given 



No{X, T) = NoiQ, k) cos[n'\Q)T - QX + 6o] , 
Ni{X, T) = iVi(g, k) cos[Q'^{Q)T - QX + 6i] . 



(46) 
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k and Q, the relative amplitudes coincides with those of the branch —k at the wavevector 
—Q, while 6{Q,k) = —6{—Q,—k). Significantly the different character of the excitations 
corresponding to the acoustic and optical branches can be seen from Fig. Elfor small wavevec- 
tors. Indeed, for the acoustic branch all functions, Nq{Q,0) ,Ni{Q,0) and N{Q,0) are of 
the same order of magnitude while the phase shift 6 ^ ivr. For the optical branches, we 
obtain No{Q, k) ^ Ni{Q, k) ^ N{Q, k) and 5 ^ vr, i.e., the concentrations of the isotropic 
and streaming electrons oscillate in antiphase, N and the space charge are almost zero. For 
finite Q, the excitations are always accompanied by a space charge wave. At A; = 1 (-1) and 
—271 <Q<0{0<Q< 27r), the phase velocity of the excitation is negative and all disturbed 
electrons concentrations are of the same order of magnitude. Otherwise, the phase velocity 
of the excitation is positive. For the developed streaming effect {Fq ^ Fcr), Ni dominates 
over A^'o. 

Thus, the accounting of space charge effects does not change considerably the basic 
properties of the electron excitations under streaming transport regime. The multi-branch 
character of the excitation spectra, the occurrence of the standing waves and other prop- 
erties can be observed at finite electron concentrations. The space charge accompanying 
the wave excitations affects mainly their damping. Space charge induced increase of damp- 
ing of electron excitations is characteristic for regimes of electron motion with considerable 
dissipation.— 

VI. DISCUSSION AND SUMMARY 

The cyclic electron motion caused by strongly inelastic processes (optical phonon emis- 
sion) is characteristic for many polar materials and associated heterostructures subjected 
to high electric field at low temperatures. The characteristic time period, Tp, and spatial 
period, Ip, of this motion depend on the applied field, F. In steady state regime, such a 
motion reveals itself via the formation of a sufficiently anisotropic distribution function and 
quasi-saturation of the drift velocity, the average energy, etc. In time (frequency) domain, 
the cyclic electron motion gives rise to the specific optical phonon transient time resonance 
at frequencies defined by the time-of-fiight, u ~ 2tt/tp. In real space, this kind of motion 
causes a spatial modulation of the electron concentration and velocity with the period Ip, 
which has been observed in structures of finite size. It is believed that the cyclic electron 
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dynamics should be more pronounced in low-dimensional structures. 

Using the BTE, we have developed an approximate approach that permits us to analyze 
coupled spatio-temporal dynamics of two-dimensional electron gas under conditions of the 
cyclic motion. As a results, we have found a novel type of excitations of the drifting electron 
gas. These wave-like excitations are periodic in time and weakly damped oscillations of the 
electrons in both real and momentum spaces. Their frequency-wavevector relations consist 
of an infinite number of continuous branches w^(g) , k = 0, ±1, ±2.... 

The specific character of the correlated motion of electrons in the wave excitations 
strongly depends on the branch number k and the wavevector q. Thus, for = , g — )■ 0, we 
have obtained oc q, the electron motion is phased, mainly, in the real space, the excitation 
is conveyed by a space charge wave. This excitation can be defined as of acoustic type. For 
k ^ and g — )■ 0, the excitations represent, mainly, strong oscillations in the momentum 
space (time-dependent redistribution between 'isotropic' and 'streaming' groups of the elec- 
trons). While the space charge is practically absent. These excitations can be defined as of 
optical type. Each of the a;'^(g)-branches once crosses the line a; = at g = 2-nk/lF- This 
implies the existence of standing waves with wavelength. A, related to the spatial period of 
the cyclic motion as A = Ip/k. For the standing waves, electron motion in real and momen- 
tum spaces is phased, but time independent. The found excitations are weakly damped. 
Their damping is caused by quasi-elastic scattering of electrons and formation of a space 
charge. The standing waves are undamped. 

Now we shall discuss the cyclic electron motion and the existence of the wave-like excita- 
tions in real polar heterostructures with two-dimensional electrons. For polar materials, the 
typical optical phonon energy is 30...100meV. Due to the exponential temperature depen- 
dence of the optical phonon absorption time (oc ^-^^o/kBT-^^ ^^n.^ process is negligible at low 
temperatures. For two examples given below, the absorption time is more than 10^ times 
larger than the emission time at T < 100 K. Other necessary and sufficient conditions for 
the realization of these effects are presented by Eqs. ([T]) and ([2]). These conditions contain 
two material parameters, the elastic scattering time in the passive region, r^, and the opti- 
cal phonon emission time, r^^. At low temperatures, the former time can be estimated by 
the use of data on the low- field mobility, /ig. The latter time can be found, for example, 
exploiting the dielectric continuum model of the optical phonons and the Frohlich electron- 
phonon interaction.— To avoid effects of electron-electron collisions on the cyclic motion. 
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the electron concentration should be less than 10^^ cm~^.— 

First, we consider AlGaN/GaN heterostructures. For two-dimensional electrons in GaN- 
channels, the optical phonon emission is estimated to be t^^ ^ 0, 01..0.02ps.^^ The high- 
est low-field mobility was measured at sub-Helium temperatures for structures with low 
density of dislocations: /io = (1...2) x lO^cm^ /Vs.^^'"^^ These results were obtained for 
sufficiently low electron concentrations: ~ (4. ..10) x 10^^ cm~^. The measurements ev- 
idence that in such structures at moderately low temperature the mobility is limited by 
acoustic phonon scattering. At T = 50. .100 K, the acoustic phonon limited mobility is 
^ (5. ..2. 5) X lO'^cm'^ /V s For numerical estimations, we choose a modest value of the 
mobility, fiQ = lO'^cm^ /Vs, then we find Tp ^ Ips. Thus the strong inequality ([T]) is met. For 
further calculations we use the following parameters: m* = 0.2mo (mo is the free electron 
mass), huo = 92meV , = 9. Then, according to Eq. @, the electric field for the onset 
of the studied effects is Fcr ~ 4.7 kV/cm. Thus, the numerical results shown in Figs. H] 
and [5] correspond to an applied field F ^ 13.7 kV/ cm. The temporal and spatial periods of 
the cyclic motion are Tp ^ 0.33 x 10^^^ s and Ip ~ 6.6 x 10~^ cm, respectively. According 
to Eq. ( 137|) . the characteristic electron concentration, for which the parameter Af = 1, is 
estimated to be N^h = 1-4 x 10^^ cm~'^. Having the scaling parameters Tp, If., N^h and using 
the results of Figs. HJ |5] one can recover the dimensional characteristics of the wave excita- 
tions. For example, the frequency of the lowest optical-like branches at g — )■ is found to be 
uj'^'^{q — )■ 0) = 1.5 X 10^'^s^^, its damping equals 7 = 6 x lO^^s^^. Undamped standing waves 
are realized for a wavevector equal to g = ±9.4 x lO^cm"^. For the same g, the frequency 
(the damping) of the acoustic and optical branches are a; ~ 2 x lO^'^s"^ (7 ~ 1.2 x lO^^s^^) 
and u; ^ 3.6 X lO^^s"^ (7 6 x lO^^s"^), respectively. 

As second example of two-dimensional systems, where the discussed effects can be ob- 
served, we consider ZnO/MgZnO heterostructures. For these strongly polar oxide het- 
erostructures, low field mobilities as high as 4 x 10^ cm? /V s were measured^^ at sub-Helium 
temperature and electron concentrations (1...4) x 10^^ cm^^. Above Helium temperature, it 
was found a phonon limited mobility temperature dependence: /xq oc T~^/'^. Particularly, 
for T ^ lOi^', it is achieved /xq ~ (2---4) x 10^ cm^/\^ s.—'^S, For numerical estimates, we 
set = lO'^cm'^ /Vs , m* = 0.3mo , huo = 72meV , hq = 8.1, = 3.7. Then, we obtain 
r^™ = 0.005...0.01ps, Tp ^ 2ps, i.e., the necessary condition of Eq. ([1]) holds. The critical 
field is Fcr = 2.9 kV/ cm. For these parameters of the ZnO/MgZnO heterostructure, results 
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presented in Figs. H] and [5] correspond to the applied field F ^ 9.7 fcV^/ cm. The periods of 
the relevant cyclic motion are Tp ~ 0.6 x 10^^^ s and Ip ~ 0.8 x 10^^ cm. The characteristic 
electron concentration N^.^ equals 0.7 x 10^^ cm~'^. With these scaling parameters and the 
results of Figs. HI O one can recover the dimensionless characteristics of the wave excitations. 

Summarizing, we have analyzed the low-temperature behavior of nonequilibrium two- 
dimensional electron gas in polar heterostructures subjected to moderately high electric 
fields. At low temperatures, when the optical phonon emission is the fastest relaxation pro- 
cess and an almost cyclic motion of individual electrons occurs, we have found the existence 
of collective wave-like excitations of the electrons. These wave-like excitations are periodic 
in time and weakly damped oscillations of the electrons in both real and momentum spaces. 
The excitation spectra are of multi-branch character, each of the spectra branches showing 
considerable spatial dispersion. There are one acoustic-type and a number of optical-type 
excitation spectra branches. Their small damping is caused by quasi-elastic scattering of 
electrons and formation of a relevant space charge. There are also waves with zero frequency 
and finite spatial periods - the standing waves. The studied excitations of the electron gas 
can be also interpreted as synchronous in time and real space manifestation of the optical 
phonon transit time resonance. Remarkably, these wave-like excitations exist in the elec- 
tron gas even in the presence of significant dissipation through optical phonon emission. 
Estimations of the parameters of the excitations for two examples of polar two-dimensional 
heterostructures, GaN/AlGaN and ZnO/MgZnO, have shown that the excitation frequen- 
cies and the standing wave periods are in THz-frequency range and sub-micrometer region, 
respectively. 
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FIG. 1: Real (a) and imaginary (b) parts of the dimensionless high frequency mobiUty of Eq. ()34p 
as functions of Q, for different Q. Results are presented for Fq/Fct = 3. Full, dashed and dotted 
lines correspond to Q = 0, 1, 2, respectively. 
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FIG. 3: Real (a) and imaginary (b) parts of the frequency as functions of the wavevector found 
from Eqs. ()40p . ()4ip for three branches of the dispersion relation at Fq/Fct = 3. Dash-dotted line 
corresponds to Vt' = VdrQ- The vertical dashed lines mark wavevectors equal to Q = ib27r. Note 
that the vertical scales in (a) and (b) are very different. 
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FIG. 4: The same as in Fig. ([3]) at = 7.4. For GaN quantum wells, this corresponds to the 
electron concentration N = 10^^ cm^^. In both panels, for comparison the results at A/" = are 
presented by the dashed lines. 
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FIG. 5: The same as in Figs. [3] and H] for the acoustic spectrum branch, k = 0, at different A^. 
For GaN quantum wehs, fuh, dash-dotted, dotted and dashed hues correspond to = 10^^, 5 x 
10^^, 2.5 X 10^1 and (in cm-^). 
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FIG. 6: Relative amplitudes and phase difference of the concentrations (see Eq. (j46p as functions 
of Q for three spectrum branches. Panels I, II, III, are for A; = 1,0, —1. In (a), the full and dashed 
lines correspond to A'^i and Nq, respectively. 
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